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Abstract 

The application of multiple-input multiple-output (MIMO) techniques to non-orthogonal multiple access 
(NOMA) systems is important to enhance the performance gains of NOMA. In this paper, a novel MIMO-NOMA 
framework for downlink and uplink transmission is proposed by applying the concept of signal alignment. By 
using stochastic geometry, closed-form analytical results are developed to facilitate the performance evaluation of 
the proposed framework for randomly deployed users and interferers. The impact of different power allocation 
strategies, such as fixed power allocation and cognitive radio inspired power allocation, on the performance of 
MIMO-NOMA is also investigated. Computer simulation results are provided to demonstrate the performance of 
the proposed framework and the accuracy of the developed analytical results. 

I. Introduction 

Non-orthogonal multiple access (NOMA) has been recognized as a spectrally efficient multiple access 
(MA) technique for the next generation of mobile networks [[I]|-[l3]|. For example, the use of NOMA has 
been recently proposed for downlink scenarios in 3rd generation partnership project long-term evolution 
(3GPP-LTE) systems, and the considering technique was termed multiuser superposition transmission 
(MUST) im. In addition, NOMA has also been identified as one of the key radio access technologies to 
increase system capacity and reduce latency in fifth generation (5G) mobile networks ifSll. If^. 

The key idea of NOMA is to exploit the power domain for multiple access, which means multiple users 
can be served concurrently at the same time, frequency, and spreading code. Instead of using water-filling 
power allocation strategies, NOMA allocates more power to the users with poorer channel conditions, 
with the aim to facilitate a balanced tradeoff between system throughput and user fairness. Initial system 
implementations of NOMA in cellular networks have demonstrated the superior spectral efficiency of 
NOMA im, O- The performance of NOMA in a network with randomly deployed single-antenna nodes 
was investigated in Q. User fairness in the context of NOMA has been addressed in |I71, where power 
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allocation was optimized under different ehannel state information (CSI) assumptions. In [[8]|, topologieal 
interferenee management has been applied for single-antenna downlink NOMA transmission. Unlike the 
above works, lIH addressed the applieation of NOMA for uplink transmission, where the problems of 
power alloeation and subcarrier alloeation were jointly optimized. The eoneept of NOMA is not limited 
to radio frequeney communieation networks, and has been reeently applied to visible light communieation 
systems in [fTOll . 

The applieation of multiple-input multiple-output (MIMO) teehnologies to NOMA is important sinee 
the use of MIMO provides additional degrees of freedom for further performanee improvement. In 
im, the multiple-input single-output seenario, where the base station had multiple antennas and users 
were equipped with a single antenna, was eonsidered. In ifT^ . a multiple-antenna base station used the 
NOMA approaeh to serve two multiple-antenna users simultaneously, where the problem of throughput 
maximization was formulated and two algorithms were proposed to solve the optimization problem. In 
many praetieal seenarios, it is preferable to serve as many users as possible in order to reduee user lateney 
and improve user fairness. Following this rationale, in [[T3ll . users were first grouped into small-size elusters, 
where NOMA was implemented for the users within one eluster and MIMO detection was used to caneel 
inter-cluster interference. Similar to mi, this method does not need CSI at the base station; however, 
unlike mi, it avoids the use of random beamforming whieh ean eause uneertainties for the quality of 
serviee (QoS) experieneed by the users. 

This paper considers a general MIMO-NOMA eommunication network where a base station is 
communieating with multiple users using the same time, frequeney, and spreading eode resources, in 
the presenee of randomly deployed interferers. The eontributions of this paper are listed as follows: 

• A general MIMO-NOMA framework whieh is applieable to both downlink and uplink transmission 
is proposed, by applying the concept of signal alignment, originally developed for multi-way relaying 
ehannels in [[T5ll and ffH . By exploiting this framework, the eonsidered multi-user MIMO-NOMA 
seenario ean be deeomposed into multiple separate single-antenna NOMA ehannels, to whieh 
eonventional NOMA protoeols ean be applied straightforwardly. 

• Since the choice of the power allocation coefficients is key to achieve a favorable throughput-fairness 
tradeoff in NOMA systems, two types of power allocation strategies are studied in this paper. The 
fixed power alloeation strategy ean realize different QoS requirements in the long term, whereas the 
eognitive radio inspired power alloeation strategy ean ensure that users’ QoS requirements are met 
instantaneously. 
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• A sophisticated approach for the user precoding/deteetion vector selection is proposed and combined 
with the signal alignment framework in order to efficiently exploit the excess degrees of freedom of 
the MIMO system. Compared to the existing MIMO-NOMA work in lfT3]l . the framework proposed 
in this paper offers two benefits. First, a larger diversity gain ean be achieved, e.g., for a scenario 
in which all nodes are equipped with M antennas, a diversity order of M is achievable, whereas 
a diversity gain of 1 is realized by the scheme in lfT3]l . Second, the proposed framework is more 
general, and also applicable to the case where the users have fewer antennas than the base station. 

• Exact expressions and asymptotie performance results are developed in order to obtain an insightful 
understanding of the proposed MIMO-NOMA framework. In particular, the outage probability is 
used as the performance criterion since it not only bounds the error probability of detection tightly, 
but also can be used to calculate the outage capacity/rate. The impact of the random locations of 
the users and the interferers is captured by applying stochastic geometry, and the diversity order is 
computed to illustrate how effieiently the degrees of freedom of the channels are used by the proposed 
framework. 

II. System Model for the Proposed MIMO-NOMA Framework 

Consider an MIMO-NOMA downlink (uplink) communieation scenario in whieh a base station is 
communicating with multiple users. The base station is equipped with M antennas and each user is 
equipped with N antennas. In this paper, we consider the scenario ^ in order to implement the 

coneept of signal alignment, an assumption more general than the one used in [[T3ll . This assumption is 
applicable to various communication scenarios, such as small cells in heterogenous networks ifTTll and 5G 
cloud radio access networks [fT^ . in which low-cost base stations are deployed with high density and it 
is reasonable to assume that the base stations have eapabilities similar to those of user handsets, sueh as 
smart phones and tablets. 

The users are assumed to be uniformly deployed in a disc, denoted by V, i.e., the cell controlled 
by the base station. The radius of the disc is r, and the base station is located at the center of V. In 
order to reduee the system load, many existing studies about NOMA have proposed to pair two users for 
the implementation of NOMA, and have demonstrated that it is ideal to pair two users whose channel 
conditions are very different |[T1, [[T9ll . Based on this insight, we assume that the disc is divided into two 
regions. The first region is a smaller dise, denoted by Vi, with radius ri (ri < r) and the base station 
located at its origin. The second region is a ring, denoted by V 2 , constructed from V by removing Vi. 
Assume that M pairs of users are selected, where user m, randomly loeated in Vi, is paired with user m'. 
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randomly located in 1^2■ Hence, the users are randomly scheduled and paired together. The use of more 
sophisticated schedulers can further improve the performance of the proposed MIMO-NOMA framework 
of course, but this is beyond the scope of this paper. 

In addition to the messages sent by the base station, the downlink NOMA users also observe signals sent 
by interference sources which are distributed in 'R? according to a homogeneous Poisson point process 
(PPP) of density A/ [|20l . The same assumption is made for the uplink case. In practice, these interferers 
can be cognitive radio transmitters, WiFi access points in LTE in the unlicensed spectrum (LTE-U), or 
transmitters from different tiers in heterogenous networks. In order to obtain tractable analytical results, it 
is assumed that the interference sources are equipped with a single antenna and use identical transmission 
powers, denoted by p/. 

Consider the use of a composite channel model with both quasi-static Rayleigh fading and large scale 


path loss. In particular, the channel matrix from the base station to user m is = 


=, where 


y/Lidm)' 

Gm denotes an x M matrix whose elements represent Rayleigh fading channel gains, dm denotes the 
distance from the base station to the user, and the resulting path loss is modelled as follows: 

, if dm > ro 


L(dm) 


Tq, otherwise 


where a denotes the path loss exponent and parameter tq avoids a singularity when the distance is small. 
It is assumed that ri > tq in order to simplify the analytical results. Eor notational simplicity, the channel 
matrix from user m to the base station is denoted by H^. Global CSI is assumed to be available at the 
users and the base station. The proposed MIMO-NOMA framework for downlink and uplink transmission 
is described in the following two subsections, respectively. 


A. Downlink MIMO-NOMA Transmission 

The base station sends the following M x 1 information-bearing vector 


s = 


+ OIm'Sm' 


( 1 ) 


where Sm is the signal intended for the m-th user, am is the power allocation coefficient, and = 1. 

The choice of the power allocation coefficients will be discussed later. 

Without loss of generality, we focus on user m, whose observation is give by 


Ym = 


Gr 


sj L(dj^ 


=Ps -f n^, 


( 2 ) 
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where P is the MxM preeoding matrix to be defined at the end of this subsection, w/^ denotes the overall 
co-channel interference received by user m, and rim denotes the noise vector. Following the classical shot 
noise model in [|2TI . the co-channel interference, w/^, can be expressed as follows: 

VPi 


A 

Wr = 

im 




liV, 


(3) 


where Im denotes an m x 1 all-one vector, and denotes the distance from user m to the y-th 
interference source. Note that small scale fading has been omitted in the interference model, since the 
effect of path loss is more dominant for interferers located far away. In addition, this simplification will 
facilitate the development of tractable analytical results. The case with p/ = 0 corresponds to the scenario 
without interference. 

User m applies a detection vector to its observation, and therefore the user’s observation can be 
re-written as follows: 


VmYm = 


G. 


yj L(ydm^ 


Ps + v^(w7^ + n^ 


(4) 


= V 


H Gm / N \ ^ H 




i^m 




Pi{aiSi + + v"(w/^ -f iim), 


interference (including inter—pair interference) + noise 

where Pm denotes the m-th column of P. 

In order to remove inter-pair interference, the following constraint has to be met: 


V^Gm 


Pi = Osxi, Vz ^ m, 


(5) 


where O^xn denotes the m x n all zero matrix. Without loss of generality, we focus on pi which needs 
to satisfy the following constraint: 


n H 


G?V2 G|{v2 


G^vm Pi — 02(m-i)xi- 


( 6 ) 


n H 


,is2(M-l)xM. 


Note that the dimension of the matrix in ®, Gfv 2 G^V 2 / G^vm G^,vm' 

Therefore, a non-zero vector pj satisfying ® does not exist. In order to ensure the existence of pi, one 
straightforward approach is to serve less user pairs, i.e., reducing the number of user pairs to (f + i). 
However, this approach will reduce the overall system throughput. 

To overcome this problem, in this paper, the concept of interference alignment is applied, which means 
the detection vectors are designed to satisfy the following constraint [[22ll . [|23]| 


^mGm — '^rn'Gm'i 


( 7 ) 
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or equivalently 


G 


H 


pH 


= 0 


Mxl- 


( 8 ) 


H 


pH 


Define as the 2N x {2N — M) matrix eontaining the {2N — M) right singular veetors of G 
corresponding to its zero singular values. Therefore, the detection vectors at the users are designed as 
follows: 


Vm 

Vm' 




(9) 


where is a {2N — M) x 1 vector to be defined later. We normalize to 2, i.e., |xp = 2, due to 
the following two reasons. First, the uplink transmission power has to be constrained as shown in the 
following subsection. Second, this facilitates the performance analysis carried out in the next section. It 


G 


H 


pH 


U X = 

^ m^m 


is straightforward to show that the choice of the detection vectors in @ satisfies 

OmxI- 

The effect of the signal alignment based design in (|7]) is the projection of the channels of the two users 
in the same pair into the same direction. Define gm — G^v^ as the effective channel vector shared by 
the two users. As a result, the number of the rows in the matrix in ® can be reduced significantly. In 
particular, the constraint for p* in ® can be rewritten as follows: 

n H 


gl 


Note that 
(fT^ exists. 

Define G = 
as follows: 


gl 

gi—1 gi+1 


gi-1 gi+1 
1 H 


gM 


Pi — O(M-l) 


xl- 


( 10 ) 


gM 


is an (M — 1) x M matrix, which means that a p* satisfying 


n H 


gl 


gM 


. A zero forcing based precoding matrix at the base station can be designed 


P = G-^D, (11) 

where D is a diagonal matrix to ensure power normalization at the base station, i.e., = 

diag{ (G-iG-ii)! 1 >" ' 1 (A)m„m dcnotcs the m-th element on the main diagonal 

of A. As a result, the transmission power at the base station can be constrained as follows: 

tr {PP^} p = tr {G"^DD^G-^} p = tr {G-^G-^D^} p = Mp, (12) 

where p denotes the transmit signal-to-noise ratio (SNR). 
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With the design in ([7]) and (fTTl) . the signal model for user m can now be written as follows: 


^mYm = 


rH 


H 


/ 

V y/L{d, 


--Pi{aiSi + ai'Si>) + v"(w/^ + n^) 


(13) 


{oimSm “f Oim'Sr 


^{L{dm))iG~^G 




+ + n« 


For notational simplicity, we define hm = 


and = 


v^rim. Therefore, the use of the signal alignment based precoding and detection matrices decomposes the 
multi-user MIMO-NOMA channels into M pairs of single-antenna NOMA channels. In particular, within 
each pair, the two users receive the following scalar observations 


2/m T Oim' ^m'') T ^mj 


(14) 


and 


(15) 


y-m' — hYa'i^CXmSm T O'm'Sm'') T T '^m'i 

where ym' and rim' are defined similar to ym and Um, respectively. Note that hm' = , \ ^ , 

and it is important to point out that hm and hm' share the same small scale fading gain with different 
distances. 

Recall that two users belonging to the same pair are selected from Vi and 1)2, respectively, which 
means that dm < dm'- Therefore, the two users from the same pair are ordered without any ambiguity, 
which simplifies the design of the power allocation coefficients, i.e., am < ttms following the NOMA 
principle. User m' decodes its message with the following signal-to-interference-plus-noise ratio (SINK) 


SINRm' = 

where the interference term is given by 


P \ hm' 1 

2^2, 

^m! 

P 1 hm' 

1 ' 1 

Vm' 

|2 + 

|v^,liv|2/m' 


T = 


(16) 


(17) 


User m carries out successive interference cancellation (SIC) by first removing the message to user m' 
with SINK, SINRmm' = i, ,2 2 , i "* 2 R'h-, \ 2 t » and then decoding its own message with SINK 

> P|ftm| a^ + \Vm\ +|V„,ljv| 1^1 


SINRm = 


|VmP + |v;^ljv|2/„ 


(18) 


which becomes the SNR if pj = 0. 
















B. Uplink MIMO-NOMA transmission 

For the NOMA uplink case, user m will send out an information bearing message Sm, and the signal 
transmitted by this user is denoted by am'VmSm- Because of the reciprocity between uplink and downlink 
channels, which was used as a downlink detection vector can be used as a precoding vector for the 
uplink scenario. Similarly P will be used as the detection matrix for the uplink case. In this paper, we 
assume that the total transmission power from one user pair is normalized as follows: 

+ awIVm'P < 2p. (19) 


The base station observes the following signal: 


M 


Ybs - 




+ 


m=l 


s/ L{dm^ sj L{^dm ') 


+ W/ + 11^5, 


where w/ is the interference term defined as follows 


A 

W/ = 


E- 




= 1m, 


( 20 ) 


( 21 ) 


je'i'/ \J L (dj.^Bs) 

dij,Bs denotes the distance between the base station and the j-th interfere^ and the noise term is defined 
similarly as in the previous section. The base station applies a detection matrix P to its observations and 
the system model at the base station can be written as follows: 


M 


PAbs = P" E 




m=l 


s/L{d„ 


s/L{d„ 


+ P "( w / + ubs)- 


As a result, the symbols from the m-th user pair can be detected based on 


V^rny BS Pri 




\ H 


\/ L{di) 


+ Pm(w/ + nBs)- 


in.terference(including inter—pair interference) + noise 

In order to avoid inter-pair interference, the following constraint needs to be met 


i^m 


Gf CTjVjSj ^ G^aj/Vj/Sj/ 


= 0, Vm 7^ i. 


( 22 ) 


Applying again the concept of signal alignment, the constraint that G^v^ = G^,Vm' is imposed on 
the precoding vectors v^. Therefore, the same design of as shown in (|9l) can be used. The total 
transmission power within one pair is given by 


pamlvml^ + pawIvm'P < pmax(a^,a^,)(|vmp + Iv^/p) < 2p. 


(23) 


Therefore, the use of the precoding vector in ® ensures that the total transmission power of one user 
pair is constrained. 
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Applying the detection matrix defined in (fTTI) . the system model for the base station to decode the 
messages from the m-th pair can be written as follows: 

yBS,m — h^O-rnSm “f Sm' + WbS + nBS (24) 

where yBS,m = WBS,m = P^w/, and nBs.m = P^t^bs- Therefore, using the proposed precoding 

and detection matrices, we can decompose the multi-user MIMO-NOMA uplink channel into M orthogonal 
single-antenna NOMA channels. Note that the variance of the noise is normalized as illustrated in the 
following: 

^{P^^BSri-lsPm] = P^Pm = (P^P)m,m = = 1- (25) 

)m,m 

The SIC strategy can be applied to decode the users’ messages, following steps similar to those used in 
the downlink scenario. 

III. Performance Analysis for Downlink MIMO-NOMA Transmission 

Two types of power allocation policies are considered in this section. One is fixed power allocation 
and the other is inspired by the cognitive ratio concept, as illustrated in the following two subsections, 
respectively. Recall that the precoding vectors and v^' are determined by as shown in @. In this 
section, a random choice of is considered first. How to find a more sophisticated choice for x^ is 
investigated in Section Illl-Ci 


A. Fixed Power Allocation 


In this case, the power allocation coefficients am and am' are constant and not related to the 
instantaneous realizations of the fading channels. We will first focus on the outage performance of user 
m'. The outage probability of user m! to decode its information is given by 


p;;, = p log 1 + 


^ Rm' I ) 


(26) 


p\hm'\‘^am + 

where P(a: < a) denotes the probability for the event x < a. The correlation between v^' and hm' 
makes the evaluation of the above outage probability very challenging. Hence, we focus on the following 
modified expression for the outage probability 


Pm/ = P log 1 -f 


p\hm'\‘^0'm + 2 + ‘251m’ 


^ Rm' 
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Since + |vmp = 2, we have < 2 and < 2. In addition, because (;^ — 

jf Y.n=i IVm'livP < iV|vm'|^- Therefore, we have 


< P. 


(27) 


for 5 > N, which means that P^' provides an upper bound on P^' if 5 > A^. Note that when 5 = 1, the 
differenee between P^' and P^' is very small as can be observed from Fig. [U i.e., a choice of 5 = 1 is 
suffieient to ensure that Pm' provides a very tight approximation to Pm' • In addition, the use of Pm' will 
be suffieient to identify the achievable diversity order of the proposed MIMO-NOMA seheme. 



Fig. 1. Comparison between P^/ and P°,, i G {m,m'}. Rm = Rm> ~ 1-5 bit per channel use (BPCU). A/ = 10 . r — 20m and 

n = 10m. ro = Im and am' = j. M = N = 2. The path loss exponent is a = 3, and the noise power is —30dBm. 

Given a random choiee of Xm, the following lemma provides an exact expression for Pm' as well as 
its high SNR approximation. 


Lemma 1. If ct^, < a^Cm', the probability Pm' = 1> where e^' = 2^'"' — 1. Otherwise the probability 
Pm' can be expressed as follows: 


~ ^ r .2 _ ^2 






<pi{x)xdx, 


(28) 


1 «/ri 


-7T\l{/3m/{x))i'y( j- 


where frn' = ^ "0 Idmfx) = 2(j)rn'SpiL (x“), and 7 (-) denotes 

the incomplete Gamma function. 

If Pi is fixed and transmit SNR p approaches infinity, the outage probability can be approximated as 
follows: 

2(/.m'(2 + L') - r“+2) 


P / ^ 

X m' 


— rt 


a + 2 


( 29 ) 













11 


where 6^’ = 2nXj6pj^. For the special case of pi = 0, P^' simplifies to 


Pm' = 1 - 




i24>rn')~ 


7 ( - + 1, - 7 f- + 1, 2(j)jn'rf 


(30) 


Proof: Please refer to Appendix A. I 

By using the high SNR approximation obtained in Lemma [T| and also the fact that both 0^' and 9„ 
are at the order of the achievable diversity gain is obtained in the following corollary. 


Corollary 1. If the diversity order achieved by the proposed MIMO-NOMA framework 

for user m' is one. 

On the other hand, user m first decodes the message for user m! before decoding its own message via 
SIC. Therefore, the outage probability at user m is given by 


P:, = P log 1 + 


P\hm 1 


P\hm\ 

1 ' 1 

|Vm 

P + |v,^l7vP/m 


^ Rm/ 


(31) 


P (log ( 1 


p\hm\‘^a'^ 


+ 


< 


Rm, log ( 1 


+ 


p \ hm 1 


P\hm 

1 ' 1 

|Vm 

P + 


Again, we focus on a modified expression for the outage probability as follows: 


> Rv 


Pm = P log 1 + 


p\hm\^oi^ + 2 + 2(5/„ 
12 ^,2 


< Rv 


(32) 


+ <Rmfiog{l + 




^ Rm' I ) 


2 + 26Im J V p\^m\‘^0(f^ + 2 + 261„ 

which is an upper bound for as explained in the proof for Lemma [2l Fig- [H demonstrates that 

Pm with a choice of <5 = 1 yields a tight upper approximation on Pm- The following lemma provides an 
exact expression for this probability as well as its high SNR approximation. 


Lemma 2. If af, < afcm', the probability Pm = 1, otherwise the probability Pm' can be expressed as 
follows: 


rri 


2 - 2 

Pm = 1 - 2 / ipi{ro)xdx - ^ I (pi{x)xdx, 


(33) 


1 JO 


1 J ro 


where fm = niax{0m, 4>m'} cmd 0m = If Pi is fixed and the transmit SNR p approaches infinity, the 

P^m 

outage probability can be approximated as follows: 


13 

-*■ m ~ 


(?)m(2 + 9m') ( a+2 , 9„a+2\ 

rlia + 2) 


( 34 ) 


where 9^' was defined in Lemma [7] 


Proof: Please refer to Appendix B. 
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B. Cognitive Radio Power Allocation 

In this section, a cognitive radio inspired power allocation strategy is studied. In particular, assume 
that user m' is viewed as a primary user in a cognitive ratio network. With orthogonal multiple access, 
the bandwidth resource occupied by user m' cannot be reused by other users, despite its poor channel 
conditions. In contrast, with NOMA, one additional user, i.e., user m, can be served simultaneously, under 
the condition that the QoS requirements of user m! can still be met. 

In particular, assume that user m! needs to achieve a target data rate of Rm', which means that the 
power allocation coefficients of NOMA need to satisfy the following constraint 


p\hjyif 1 



|Vm' 

|2 + 

|v^,l7vP/m' 


^ Cm') 


which leads to the following choice for a. 

“m = max ( 0, 


P| ^m' 

P - em'(|Vm'P + h 


to 

(1 + ^m')p\hm‘ 

,|2 


(35) 


(36) 


It is straightforward to show that 




is always less than one. 


An outage at user m' means here that all power is allocated to user m', but outage still occurs. As 
a result, the outage probability of user m' is exactly the same as that in conventional orthogonal MA 
systems. Therefore, in this section, we only focus on the outage probability of user m which can be 
expressed as follows: 


=P (|/i^p < max {0m'(|v^p + Iv^ljvp/^), + 


I 


■w}). 


(37) 


if Q(^, > otherwise outage always occurs. It can be verified that a^, < a^Cm' is equivalent to 

am = 0, in the context of cognitive radio power allocation. 

Analyzing this outage probability is very difficult due to the following two reasons. First, hm and 
Vm are correlated, and second, the users experience different but correlated co-channel interference, i.e., 
Im 7^ Ira'- Therefore, in this subsection, we only focus on the case without co-channel interference, i.e., 
p/ = 0. In particular, we focus on the following outage probability 

Pm =P (|/imP < 2max{(^m/,^^}) , (38) 


where chm = 


pa'C-palntm' 


and 


Qfm = max ( 0, 


p|/im'| 26m' 


(39) 


(1 T Cm')p|^m'P, 

Similarly to the case with fixed power allocation, the outage probability Pm tightly bounds Pm- The 
following lemma provides the expression for the outage probability Pm- 
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Lemma 3. When pi = 0, the outage probability can be expressed as follows: 


Pm = 1 - Ti 


26m' 

P 


T. 


26)71 (1 “ 1 “ 6 ) 7 )/ ) 


where 


Ti(|/) = 




— rr 


P 


7 - + 1,1/r" 

a 


7 - + l,l/r“ 

« 




Tofz') = 


r^e-^^o 






“-0^) + 


a 


7 - + l, 2 :r“ - 7 - + l, 2 :r, 


a 


At high SNR, the outage probability can be approximated as follows: 


46)77' 


_ ^«+2 


r) + 


2r2+“ 


6m(1 + 6,n') ^ 46m(1 + ^m') (^a+2 ^a+2 

/^ . N o V'l 'l 


pr? 


p(2 + a)rl 


(40) 


(41) 


) . (42) 


p{2 + a){r‘^ — rf) 

Proof: Please refer to Appendix C. ■ 

Remark 1: By using the above lemma, it is straightforward to show that a diversity gain of one is still 
aehievable at user m (i.e., there is no error floor), and it is important to point out that this is achieved 
when user m' experiences the same outage performance as if it solely uses the channel. Therefore, by 
using the proposed cognitive radio NOMA, one additional user, user m, is introduced into the system to 
share the spectrum with the primary user, user m', without causing any performance degradation at user 


m'. 


Remark 2: For the above cognitive radio NOMA scheme, it was assumed that the message for user m' 
is decoded first at both receivers. Nevertheless, different SIC decoding strategies can be used, and their 
impact can be obtained in a straightforward manner from the analysis in the next section, where more 
complicated uplink transmission schemes are studied. It is worth pointing out that in (1^ is always 
smaller than 4, for Rm' > 1. For example, when = 0, the inequality — 4 < 0 holds obviously. 
When > 0, 


2 1 p\hra' 


p + h 

1 p\hm'\ 

(1 6)77') 26)77' ( 1^)77 1 

- _L 1 |2 r ^ 

«m 2 

(1 “h ^m'^p\h7n‘ 

'|2 2 

2(1 ^m'')p\hm^ 

,2 


if Rm' > 1. 


( 43 ) 
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Algorithm 1 The selection of the detection vectors v^, and 


1 : for i = 1 to {2N — M) do 


2: Set ^rri.A 


3: 

4: 


Oix(i-l) 1 Oix(Ar-i) 


H 


, Vm e {1, • • • , M}. 


Choose the detection vector as 


n H 


■ V"^/ • 

m,^ m',2 


= UmXm,i and determine vector 


Construct the effective small scale fading matrix, denoted by Gj, by using gm,i, i.e., G* = 

H 


gl,i 


gM,) 


5: Find the effective small scale fading gain for each user pair, 'jrn,i = 

6: Find the smallest fading gain, 7niin,i = minjyi^j, • • • , 7 m, i}- 

7: end for 

8: Find the index i which maximizes the smallest fading gain, i* = 




arg max 7min,i ■ 

,2N-M} 


C. Selection of the User Detection Vectors 

Previously, a random choice of and v^' has been used and analyzed. In the case of 2N — M > 1, 

there is more than one possible choice based on the null space, Um, defined in (|9l). In this section, we 

study how to utilize these additional degrees of freedom and analyze their impact on the outage probability. 

Finding the optimal choice for and v^' is challenging, since the choice of the detection vectors for 

one user pair has an impact on those of the other user pairs. For example, the choice of and v^' will 

affect the m-th column of the effective fading matrix G. Recall that the data rates of the users from the 

i-th pair is a function of • Therefore, the detection vector chosen by the m-th user pair will 

also affect the data rates of the users in the i-th pair, m i. 

In order to avoid this tangled effect, a simple algorithm for detection vector selection is proposed in 

Table [B The following lemma shows the diversity gain achieved by the proposed selection algorithm. 
Lemma 4. Consider the use of a fixed set of power allocation coefficients. Ifcf^, < the probability 

Pm' = 1, otherwise the use of the algorithm proposed in Table\I\ ensures that a diversity gain of {2N — M) 

is achieved. 

Proof: Please refer to Appendix D. ■ 

As can be seen from Lemma IH the use of the proposed selection algorithm can increase the diversity 
gain from 1 to {2N — M), which is a significant improvment compared to the scheme in lfT3ll . Consider a 
scenario with N = M as an example. The proposed scheme can achieve a diversity gain of M, whereas 
the one in [fT3ll can only achieve a diversity gain of 1, for an unordered user. Note, however, that the 
scheme in lfT3]l does not require CSI at the transmitter. 
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IV. Performance Analysis of MIMO-NOMA Uplink Transmission 

Because of the symmetry between the uplink and downlink system models of Section [III in this section, 
we only focus on the difference between two scenarios. One important observation for uplink NOMA is 
that the sum rate is always the same, no matter which decoding order is used. Therefore, in this section, 
we first analyze the outage probability with respect to the sum rate for a fixed power allocation. The use 
of a randomly selected is considered in order to obtain tractable analytical results. 


A. Fixed Power Allocation 


Recall that, if the message from user m is decoded first, the base station can correctly decode the 
message with rate 


R 


m^BS. 


/ = log 1 + 


P I I ^ ^BS,m T 1/ 


where the interference power is given by 


I 


BS,m 


E 


Pi \ Pm 


LM 


(44) 


(45) 


After subtracting the message from user m, the base station can decode the message from user m! 
correctly with the following rate 


(47) 


Rm',BS,I = log 1 + 

V lBS,m + J- 

Therefore, the sum rate achieved by NOMA in the m-th sub-channel is given by 


= R 


m,BS,I + Rm',BS, 


3,1 = log 1 


7sS,m + 1 

It is straightforward to verify that the exactly same sum rate is achieved if the message from user m' is 
decoded first. Therefore, the outage probability for the sum rate can be expressed as follows: 


Pg — P (i?s < + Rji 


(48) 


Note that the term for the interference power contains which makes the calculation very difficult. 

Since |p^1mP < we focus on the following modified expression of the outage probability 


P^ = P log 1 + 


+ 1 


1 ^ p _1_ P 


(49) 


where I^= Y, 


PI 




. Similarly to the downlink case, P^ provides an upper bound on P^ for 6 > M. In 


the simulation section, we will demonstrate that P^ with a choice of 5 = 1 provides a tight approximation 
to P^. 
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Define the small scale fading gain as x = 
expressed as follows 


1 




The sum rate outage probability can be 


p p, ,1 pL., f(«„ + i) 


Sim + 1 




(50) 


L{dm) ' L(d^,) 

where e = 2^’"+^’"' — 1. Following the same steps as in the proof of Lemma [H the above probability can 
be expressed as follows: 


p =1_!_ r r 

- rl) I, io 




xdxydy, 


(51) 


where C{dm,dm') = 

L(dm) ) 

In order to obtain some insights regarding the above probability, we again consider the case that p tends 
to infinity and pi is fixed. Since both dm and dm’ are bounded, ({dm, dm') approaches zero at high SNR. 
Therefore the above probability can be approximated as follows: 


1 - 

1 - 



e-a-,y)e~''^dpisax,y ))« 


3 ro ^^)xdxydy. 


xdxydy 


(52) 


With some algebraic manipulations, the above probability can be simplified as follows: 

TrXjSotpj 


T3 

1 S ^ 



1 - ({x,y) 


({x, y)xdxydy. 


To 


+1 xdxydy 


(53) 


Therefore, the outage probability can be approximated as follows: 

4{e {’’^ + l) 1 


s 


(54) 


where ^ = T f, 


prf(r^ — rf) p' 

r dxdy is a constant and not related to the SNR. Hence, a diversity gain of 1 is 


ri Jo 2 

I T7i‘ 


L(x) L{y) 

achievable for the sum rate. 


B. Cognitive Radio Power Allocation 

The design of cognitive radio NOMA for uplink transmission is more complicated, as explained in the 
following. To simplify the illustration, we omit the interference term in this section, i.e., pi = 0. For 
downlink transmission, < 4 was sufficient to decide the SIC decoding order. However, there are more 
uncertainties in the uplink case, since is not necessarily larger than al^\hm\^ even if > 4. 
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Therefore, the base station ean apply two types of deeoding strategies, i.e., it may deeode the message 
from user m' first, or that of user m first. These strategies will yield different tradeoffs between the outage 
performanee of the two users, as explained in the following subseetions, respeetively. 

1) Case I: When the message from user m! is decoded first, in order to guarantee the QoS at user m', 
we impose the following power constraint for the power allocation coefficients 


log(l+ 


which leads to the following choice for 


“m' = min 1, 


p\hm\^al^ + 1 / ^ ’ 


+ pe 

m' I 


(55) 


(56) 


+ em'P\hm\‘^ 

Following the same steps as in the proof of Lemma [2l the outage probability P^, can be evaluated 
as follows: 

+ pe m' I hm I 




p\hm\ 

^m/ 


> 1 


(57) 


= p 


p 


= 1 - Ti 


P 


and following the same steps as in the proof of Lemma [H the outage probability can be evaluated 

as follows: 

T P^m'\hm\ 


=P 


^m' T P^m' I h m I 


> 1 


(58) 


p\hm'\'^ T ^m'P\hm\'^ 




< l,log 1 + p\h 


=P +P X> 


=1 - Ti 


P 

Cm' (f T Cm) 
P 


12 P|(im'| Cm' ^ 

P|/im'P + Cm'P|(imP/ 

, X < L[dm,) - h L[dm') - (1 + Cm) 

P P P 


^m' L{^dm'') 


(59) 


T, ( y 1 . 


2) Case II: When the message from user m is decoded first, in order to guarantee the QoS at user m', 
we impose the following power constraint for the power allocation coefficients 

log (l + p\hyn'\^al^) > Rm', (60) 

which leads to the following choice for ctm' 


«m' = min 1, 




p\K 


(61) 


With this choice, we can ensure that the outage probabilities of both users are identical, i.e., P^^ 


m,BS 


Pm' explained in the following. The outage events that occur at user m' can be divided into the 

following three events 
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• El'. All the power is allocated to user m', i.e., a^' = 1, but the user is still in outage. The NOMA 
system is degraded to a scenario in which only user m' is served. 

. E 2 '. When < 1, outage occurs at user m, and SIC is stopped. 

• E^i'. When < 1, no outage occurs at user m, but outage occurs at user m'. 

It is straightforward to show that E^ will not happen, i.e., P(.E 3 ) = 0. Therefore Prn',BS = P(-^i) +P(-^ 2 )- 
On the other hand, there are only two outage events for decoding the message from user m, which are Ei 
and E 2 , respectively. Therefore, the outage probabilities of the two users are the same, P^ ss = Pm' bs- 
Therefore, we only need to study the outage probability for the message from user m. With the choice 
shown in (IMl) . the outage probability can be rewritten as follows: 


pii _ p 
^m,BS — ^ 




p\hr 


> 1 +P 


p\h 


< 1 




p\h 


, |2_ 'LmL _I- 1 

■m'\ p\h^,\2 + ^ 




(62) 


Therefore, the outage probability can be expressed as follows: 




.,BS 


= P X < 


P(*^m')Cn 

P 


P 


+ 1) L[d 

m) 

OC 


p 


p 


(63) 


By applying the same steps as in the proof of Lemma[3]for finding P(i?i) and P(T^ 3 ), the outage probability 
can be obtained as follows: 

em(l + Cm') 


pll _ pll _ 1 
^m',BS — ^m,BS ~ 


Ti ( ^ 

P 




P 


(64) 


Remark 3: The two considered cases strike different tradeoffs between the outage performance of the 
two users. Case I can ensure that the QoS at user m' is strictly met, and therefore user m' will experience 
a lower outage probability in Case I, which can be confirmed by the fact that Pm',B 5 < Pm',B 5 ’ due to 
T 2 j < X On the other hand. Case II does not require that the message of user m! arrives at 

the base station with a stronger signal strength since the base station will decode the message from user 
m first. This is important to avoid the problem of using too much power for compensating the huge path 
loss of the channel of user m!. As a result, more power is allcoated to user m compared to Case I, and 
hence, user m experiences better outage performance in Case II, i.e., P^ ^5 > P^ bs- This can be shown 
by comparing (l59l) with (1^ and by considering 


T(dm)Cm T £m' (^m T 1) T(dm') ^ Li{,dm'^^m' T (^m' T 1) T((im)- 


(65) 


V. Numerical Studies 

In this section, the performance of the proposed NOMA framework is investigated by using computer 
simulations. The performance of three benchmark schemes, termed MIMO-OMA without precoding. 
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MIMO-OMA with precoding, and MIMO-NOMA without precoding, is shown in Fig. |2l in order to better 
illustrate the performance gain of the proposed framework. The design for the two schemes without 
precoding can be found in [[T3]l . The MIMO-OMA scheme with precoding serves M users during each 
orthogonal channel use, e.g., one time slot, whereas 2M users are served simultaneously by the proposed 
scheme. For MIMO-OMA with precoding, the design of the detection vectors was obtained by following 
the algorithm proposed in Table [U where the users will carry out antenna selection in each iteration. The 
framework proposed in this paper is termed SA-MIMO-NOMA. The path loss exponent is set as a = 3. 
The size of Vi and V 2 is determined by r = 20m, and ri = 10m. The parameter for the bounded path 
loss model is set as Tq = 1. 




(a) Outage Sum rate (b) Outage Probabilities 

Fig. 2. Performance comparison with the three benchmark schemes for downlink transmission. Rm = 5 BPCU and Rm' = 0.5 BPCU. 
r = 20m and ri = 10m. M = N = 3. ro = Im. a„,/ = |. The path loss exponent is a = 3. The noise power is —30dBm and the 
interference power is pi = 0. 


Since the benchmark schemes were proposed for the interference-free scenario. Fig. [2] shows the 


performance comparison of the four schemes for p/ = 0. In Fig. |2(a)[ the downlink outage sum rate, 
defined as Pm), is shown as a function of transmission power, and the corresponding 


outage probabilities are studied in Fig. |2(b)[ As can be seen from the figures, the two NOMA schemes 
can achieve larger outage sum rates compared to the two OMA schemes, which demonstrates the superior 
spectral efficiency of NOMA. In Fig. |2(b)[ the two schemes with precoding can achieve better outage 
performance than the two schemes without precoding, due to the efficient use of the degrees of freedom at 
the base station. Comparing SA-MIMO-NOMA with the MIMO-NOMA scheme proposed in lfT3]l . one can 
observe that their outage sum rate performances are similar, but SA-MIMO-NOMA can offer much better 
reception reliability, particularly with high transmission power. In terms of individual outage probability. 
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SA-MIMO-NOMA can ensure a lower outage probability at user m, i.e., a smaller P^, compared to the 
MIMO-OMA scheme with precoding, but results in performanee degradation for the outage probability at 
user m', i.e., an inerease of P^'- This is eonsistent with the finding in [|T^ which shows that the NOMA 
user with poorer ehannel conditions will suffer some performanee loss due to the co-ehannel interference 
from its partner. 




(a) Exact expressions (Rm' = Rm = 1 BPCU) (b) Approximation (pi = — lOdBm) 

Fig. 3. Outage probabilities V^n' and Vm for downlink transmission. Xi — 10~^, <5 = 1, r = 20m, n = 10m, M = N = 2, ro = Im, 
and a„i/ = The path loss exponent is a = 3 and the noise power is —30dBm. The analytical results are based on Lemmas [T] and 



Fig. 4. Outage probability for cognitive radio downlink transmission, r — 20m, ri = 10m, ro = Im, S = 1, pi = 0, and M = N = 2. 
The noise power is —30dBm. The analytical results and the approximations are based on Lemma 


In Fig.[3l the accuracy of the analytieal results developed in Lemmas [U and [21 for downlink transmission 
is verified. As ean be seen from Fig. |3(a)[ the exaet expression developed in Lemma [H perfeetly matches 
the computer simulations, and the asymptotic results developed in Lemma [H are also aeeurate at high 
SNR, as shown in Fig. |3(b)[ The aeeuracy of Lemma [2] can be eonfirmed similarly. Note that error floors 
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appear when inereasing pj in Fig. |3(a)[ which is expected due to the strong co-channel interference caused 
by the randomly deployed interferers. 

In Fig. m the performance of the cognitive radio power allocation scheme proposed in Section ITTT-RI 
is studied. In particular, given the target data rate at user m!, the power allocation coefficients can be 
calculated opportunistically according to (l3^ . As can be seen from the figure, the probability for this 
NOMA system to support the secondary user, i.e., user m, with a target data rate of Rm approaches one 
at high SNR. Note that with OMA, user m cannot be admitted into the channel occupied by user m!, and 
with cognitive radio NOMA, one additional user, user m, can be served without degrading the outage 
performance of the primary user, i.e., user m'. 

In Fig. [5l the impact of the number of user antennas on the outage probability is studied. As can be 
seen from the figure, by increasing the number of the user antennas, the outage probability is decreased, 
since the dimension of the null space, Um, defined in (|9l), is increased and there are more possible choices 
for the detection vectors. Furthermore, the slope of the outage curves is also increased, which indicates 
an increase of the achieved diversity order and hence confirms the findings of Lemma 01 



Fig. 5. Impact of the number of the user antennas on the downlink outage probabilities Pm and Pm' - Rm = 4 BPCU, Rm' = 1-9 BPCU, 
Om' = |, pr = 0, r = 20m, ri = 10m, and rg = Im. The noise power is —30dBm. 

The performance of the proposed NOMA framework for uplink transmission is demonstrated in Figs. 
[6] and |71 In particular, in Fig. the outage probability for the sum rate is investigated, and in Fig. |7] 
the performance of the proposed cognitive radio uplink schemes is studied. As can be observed from 
both figures, the developed analytical results perfectly match the computer simulation results, which 
demonstrates the accuracy of the developed analytical framework. It is worth pointing out that the modified 
probability P* with 5 = 1 provides an accurate approximation for P^. An interesting observation from 
Fig. |7] is that Cases I and II offer different performance advantages. In terms of P^s Case I can offer 
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Fig. 6. Outage probabilities for the uplink sum rate and Ps- r = 20m, ri = 10m, tq = Im, 5 = 1, A/ = 10“^, pi = lOdB, and 
M = N = 2. The noise power is —30dBm, and the interference power is pi = — lOdBm. The analytical results and the approximations are 
based on (ED and El, respectively. 



Fig. 7. Uplink outage probability for user m with the cognitive radio constraint, r = 4m, ri = 2m, ro = Im, pi — 0, and the noise 
power is —30dBm. The analytical results and the approximations are based on ( 158b and El, respectively. 

a lower outage probability compared to Case II, however it results in a loss in outage performance for 
user m. In practice, if the QoS requirement at user m' is strict. Case I should be used, since the outage 
probability realized by Case I is exactly the same as when the entire bandwidth is solely occupied by 
user m'. Otherwise, the use of Case II is more preferable since the outage performance for user m can 
be improved and the system will not spend exceedingly high powers to compensate the user with poorer 
channel conditions. One can also observe that, for Case I with Rm' = Rm, the outage performance for 
user m is worse than that of user m!, although user m is closer to the base station. The reason for this is 
because in Case I, the power is allocated to user m' first, and user m is served only if there is any power 
left. Therefore, the outage probability of user m will be at least the same as that of user m!, as discussed 
in Section |IVl 
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VI. Conclusions 

In this paper, we have proposed a signal alignment based framework whieh is applieable to both MIMO- 
NOMA downlink and uplink transmission. By applying tools from stoehastie geometry, the impaet of the 
random loeations of the users and interferers has been eaptured, and elosed-form expressions for the 
outage probability aehieved by the proposed framework have been developed to faeilitate performance 
evaluation. In addition to fixed power allocation, a more opportunistic power allocation strategy inspired 
by cognitive ratio networks has also been investigated. Compared to the existing MIMO-NOMA work, the 
proposed framework is not only more general, i.e., applicable to both uplink and downlink transmissions, 
but also offers a significant performance gain in terms of reception reliability. In this paper, it has been 
assumed that global CSI is available, which may introduce a significant training overhead in practice. An 
important future direction is to study how MIMO-NOMA transmission can be realized with limited CSI 
feedback. 


Appendix A 
Proof for Lemma □ 

First, we rewrite the considered probability P^' as follows: 

-v2 


pa 


Prr>' = P 




_ 

L{d^,){G-^G-^)r. 




( 66 ) 


+ 2 + ‘251m' 

In order to calculate P^s the density functions for the three parameters, dm', Im' and — 

have to be found. Recall that the factor ^ 


1 


(G-iQ 




can be written as follows m-. 

~ Sm {'^M ©m) gmj 


(67) 


/m,m 


where Qm = Gm(G^Gm)~^G^ and G^ is obtained from G by removing its m-th row. If gm is complex 

Gaussian distributed, the density function of (q^iq-h-) - will be exponentially distributed. This can be 

shown as follows. First, note that the projection matrix (Im — ©m) is an idempotent matrix and has 
eigenvalues which are either zero or one. Second, recall that each row of G is generated from an M x 2N 


complex Gaussian matrix 




gm 2 


, I.e., 


'^m '^m' 


m m' 


H 




Id m, X!rr7 


( 68 ) 


Hence, provided that is a randomly generated and normalized vector, the application of Proposition 1 
in [l2^ yields the following 


gri 


CN(0,Im), 


(69) 
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i.e., gm is still an M x 1 complex Gaussian (CN) vector. Therefore, -is indeed exponentially 

^ )m,7n 

distributed, and the outage probability can be expressed as follows: 




Hd'm' ) Hd-m' )^m' 

'-V-' 

Qi 


(70) 


which is conditioned on Otherwise, Pm' is always one. 

Since the homogenous PPP \h/ is stationary, the statistics of the interference seen by user m' is the 
same as that seen by any other receiver, according to Slivnyak’s theorem Il25]l . Therefore, Im' can be 
equivalently evaluated by focusing on the interference reception seen at a node located at the origin, 
denoted by Jq = \ » where denotes the distance between the origin and the j-th interference 

ie'E'/V _ _ 

source. As a result, the expectation of Qi with respect to Im' can be expressed as follows: GOll . Il26ll 

£i„, {Qi} = £i„, < e I = exp (1 - dpj , (71) 

where p denotes the coordinate of the interference source, and d denotes the distance. Note that distance d 
is determined by the node location p. After changing to polar coordinates, the factor Sj , can be calculated 


as follows: 

{Qi} = exp ( -TiXirl ( 1 - e“ 


exp ( —27rA7 




xdx 


= exp -T^XlPm'l 


1 A 


a r. 


0 


(72) 


where Pm'{dm') is denoted by Am' for notational simplicity. Therefore, the outage probability can be 
expressed as follows: 




(73) 


Recall that user m! is uniformly distributed in the ring V 2 . Therefore, the above expectation with respect 


to dm' can be calculated as follows: 


Pm' = l- 

dpGT>2 


dp 


9 9 ’ 


(74) 


— vrrj 

where distance dm' is determined by the user location p. Changing again to polar coordinates, this 
probability can be expressed as follows: 

2 


Pm' = 1 2 2 






(75) 


1 J ri 


Hence, the first part of the lemma is proved. 
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In the case that p approaches infinity and pj is fixed, it is easy to verify that (pm', as well as jSm', go 
to zero. Hence, the incomplete Gamma function in (17^ can be approximated as follows: 


7 


1 Pm' \ _ 

rv ’ T? ; Z ^ 


a Tn 




n\ 


(^ + ^) 


a 


^ ^ n=0 

Therefore, the factor 8i , can be approximated as follows: 


Pn 


(76) 


{Qi} ~ exp -nXiP^,a 


Pn 


— g Tji' m ^ 


(77) 


where 6m' = 27iXj6pm'Pi-^- Using this approximation the outage probability can be simplified at high 


SNR as follows: 


Vm' ~ 1-^ / e 


2 _ .^2 


tyt 


(78) 


l J ri 


2 _ ^2 




(1 — {2pm' + 6m')x°) Xdx = 


I J ri 


2{2Pm' + dm') 

r'^ — r? a + 2 


For the special cause without co-channel interfere, i.e., pj = 0, the probability in (1751) can be simplified 
as follows: 


Pm' = 1- f e ‘^^'"'^‘^xdx = 1 


2 _ .^2 




1 «/ri 


/p2, _ /jo2 / 




(79) 


1 6 ( -r°‘ 2 -rf 2^ (2^m') 

= 1-^- 2 (e r -e ^ 

1 




?) 


2 _ 2 \1 + ^,‘2pm'r°‘] - 7 (z+ ^,‘2pm'ri 




1 


a 


a 


and the lemma is proved. 


Appendix B 
Proof for Lemma [2] 

When > op^em', the outage probability can be written as follows: 

Pm = P {\hmP < 2pm'{^ + Sim)) + P (|/im|^ < 2pm{^ + dim): \hmP > 20m'(l + dim)) (80) 
= P {\hmP < 2m&-\{pm:Pm'}{^ + dim)) : 

The reason why P^ is an upper bound on Vm for 5 can be explained as follows. Recall that the orig¬ 
inal outage probability P^ can be expressed as P^ = P [\hmP < max{0m, 0m'}(|vmP + |v;^lArp/m))- 
Since |v^l 7 vp < and |vmp < 2, we have Vm < Pm if 5 > iV. It is worth pointing out that a 

choice of 5 = 1 is sufficient to yield a tight approximation on P^, as shown in Fig. [T] 

Recall that hm = , \ _ =■ Comparing hm to hm', we find that the only difference between 

the two is the distance dm which is less than ri. In addition, the statistics of Im can be studied by using Iq 
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as explained in the proof of Lemma [IJ Therefore, following steps similar to those in the proof of Lemma 
[H the outage probability ean be expressed as follows: 

P _ O \ X _ 6 2077],Z/(rf77T,) 

7Tl | 

[ Q2 

It is straightforward to show that the expeetation of <52 can be obtained in the same way as that of Qi, by 
replaeing cpm' with (pm- In addition, recall that user m is uniformly distributed in the disc Vi. Therefore, 
the outage probability can be calculated as follows: 



Pm = l- [ (82) 

Jpev, Trrf 

where distance dm is again determined by the user location p. Resorting to polar coordinates, the outage 
probability can be expressed as follows: 

2 - 2 

Pm = 1- 5 - / ipi{L{x)xdx - 5 - / e~‘^‘^"'^°‘(pi(L(x))xdx. (83) 

n Jo n Jro 

If p approaches infinity and pj is fixed, both Pm and pm go to zero. With this approximation, the incomplete 
Gamma function in (1721) can be approximated as {Qi} ~ where 9m = 2Ti\ipmpif^- Hence, 

the outage probability can be simplified at high SNR as follows: 


p ~ 1 — 

2 r° 
~ / ^ 

> '^4>mrQ ^ 

a ^0 'Tfl^ 

^ kaj 


rf Jo 



^ 1 — 

4(i- 

2 pmrQ 

- r^Om^ - 


rf \ 

^( 24 ,.. 

rf{ 

m ^m) 

a + 2) 

(ar“+2 

+ 2r“+2) ; 


fri 


n Jr, 

2 n 


r 


1 Jr, 


1 — {2pm + 9m)x‘^ ) xdx 


(84) 


and the lemma is proved. 


Appendix C 
Proof for Lemma [3] 

There are three types of outage events at user m, as illustrated in the following: 

• am = 0, i.e., all the power is consumed by user m' and no power is allocated to user m. This event 
is denoted by Ei. 

• When > 0, user m cannot decode the message to user m' . This event is denoted by E 2 . 

• When > 0, user m can decode the message to user m', but fails to decode its own message. This 

event is denoted by E^^. 
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The probability of Ei can be expressed as follows: 

P(Ei) = P {p\hm'\^ - 2em' < o) . 


(85) 


This probability can be straightforwardly obtained from the proof of Lemma [T| by replacing cpm' with 


A 2e, 


Therefore, P(-E'i) can be expressed as follows: 


P(,Bi) = 1 - Ti 




When am > 0, P(i? 2 ) = 0, since 

+2 ^ + 2)) 

=P {^p\hmf < a^p\hmf{^ + Cm') + 26^') 

=P {p\hm\^\hm'? < p\hm'?\hm? “ 2|/i^pe^/ + 2\hm'?em') = P {\hm\^ < \hm'?) = 0- 


( 86 ) 


(87) 


The probability for event i ?3 can be calculated as follows: 


P(E 3 ) = P (log (^1 + < Rm,Oim > 0 j = P H/l, 


12 P|(^m'| 2€.m' I, |2 ^ 26^' 

^ ) I ^m' \ ^ 

P 


2 (1 T €m ') I hm' I 


( 88 ) 


An important observation is that both channel gains hm and hm' share the same small scale fading. 
Defining x = -, the outage probability can be expressed as follows: 

2er; 


P(g3) = p| .7, " 


> 


L{dm) (1 + ^m’) ,) L{dm') P 

T-, I L{dm'^ 2,(im'L{^dm'^ 26m(l T 

= P I - < X < - 1 - 

p 


(89) 


P P 

The above probability can be calculated as follows 


P(i?3) = 


'dpr, 



dpmdpm' 1 


(90) 


Pm '^ T^2 


PmeX>l,p„,/GX>2 


where pm denotes the location of user m. Since the users are uniformly distributed, the above probability 
can be expressed as follows: 


P(B3) = 




fP 2e , 

_ m 

e py°‘ ydy — 


1) Jri 


r?(r^ — r 


1 ) Jo 


2em(l + e„./) 2e^iL{x) 

e pv°‘ ydy e p xdx 
J ri 


(91) 


= Ti 


2eri 


P 


-Ti 


26,; 


Tc 


2e 

ra (1 + ) 

p 
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Combining (l86l) . ([87]) . and (ITO . the first part of the lemma ean be proved. To obtain the high SNR 
approximation, we have 


Ti(|/) ^ 1 + ^-2 + 


y~ 




= 1 - 


2 ?/ 


(2 + a){r‘^ — rf) 

when y approaehes zero, and 






(92) 


T2{z) ^ 1 - 


„2+a. 
^0 




oi.-\-2 _ q :+2 

^ / I 


) + 




(i + l)rf 


' ct \ — 
ZTi j “ 


[zr^ 


>4+1 


= 1 - 




2z 


(2 + a)rl 


(r“+2 _ 


(93) 


when z approaehes zero. By substituting the above approximations into (l40l) . the lemma is proved. 


Appendix D 
Proof for Lemma H] 

We foeus on the outage performanee of user m! first. Given the deteetion veetor \rn,i* chosen from 
Table [B the outage probability ean be upper bounded as follows: 


/ 


Pm',i* < P 


pa^ , 




\ 




y Lid^,){G-^G-,^)m.r, 

= P (7m,i* < 2(l)m'L{dm'){l + dim')) < P (7mm,i* < 2(t)m'L{dm'){l + dim')) ■ (94) 

According to the algorithm proposed in Table [T1 


+ 2 + 2dl„ 




/ 


7min,i* — 1113^{7min,l; ‘ ‘ ‘ i 7min,2 A'—m}- 


(95) 


s2N-M 


Therefore, the outage probability can be bounded as follows: 

Pm',i* ^ (P (7min,i 7 2(f)m'L{dm'){^ T dim'))) 

where the inequality follows from the fact that 7min,i and 7min,i are independent, since and gmj are 
independent (Proposition 1 in lf23]l l. The above outage probability can be further bounded as follows: 


Fm',i* < (37P (7m,i < 2(pm'L{dm'){l + dim'))) 


2N-M 


(96) 


Following the same steps as in the proof of Lemma [B the upper bound on the outage probability can 
be calculated as follows: 


Pm',i* < |(1 - 


2N-M 


< M' 


2N-M 


2N-M 

i=0 


2N-M 




( 97 ) 
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which is conditioned on 


After the expeetation with respeet to the outage probability ean be bounded as follows: 


Pm',** < M' 


2N-M 


2N-M 

E 

*=o 


2N-M 


i-iySd , < 




(98) 


For the ease of p approaehing infinity and a fixed pj, the upper bound on the outage probability ean 
be approximated as follows: 


2N-M 




< M' 


2N-M 


i=0 

2N-M 

E 

*=o 


2N -M' 


2N -M 


{-irSd, { 




-7rAj(j/3^/)Sa( 


(99) 


Using polar coordinates, the upper bound ean be calculated as follows: 

2N-M 


Pm',*‘ < Y 


*=0 


2N-M 


J.2 _ J.2 




■xdx (100) 


1 i=0 


2M2 


_ ^*>2 


E 

i=0 


■i)‘E 

j=0 


2N-M\. ,, ^ (-!)'■ (*^m' + 


j! 


ia + 2 


( 101 ) 


By exchanging the two sums in the above equation, the upper bound ean be rewritten as follows: 


P m/ i* ^ 


2M' 


2N-M °° 


m', 1 * ^2 ^2 


E 


(-l)'(9m' + (7“+2 - rf+■') /2]V - M 


1 


2M' 


j=0 
2N-M °o 


2 ! 


ia + 2 


E 

*=o 




( 102 ) 


2 _ 




1 „• 


E 


(-1)J(6>^. + 2(Pm'y (2N - M' 


j=2N-M 




where the last step follows from the following fact 


2N-M 

E 

i=0 


2N-M 


ja + 2 


(-l)V = 0 


E 

i=Q 


(103) 


for 0 < j < (2A^ — M — 1) GTll . Furthermore, note that both and 9m' approaeh zero for the considered 
seenario, and Ym=o~^ ^_^y^ 2 N-M _ — M)\. Therefore, the upper bound on the 

outage probability can be approximated as follows: 

t2N-m ( T\2N-Mra I oA ^2Af-M | 


Pm',i* E 


+2(j)m'y 


<^2 _ ,^2 


{2N-M)\ {2N-M)a + 2 

2[M{9m' + 2cl>m')?^-^ ( r(2^-^)“+2 _ ^(2A^-M)a+2> 


- M)\ 


( 104 ) 


(r2-r2)((2A^-M)a + 2) 


P 


2N-M ■ 


The result for user m ean be proved using steps similar to the ones above. 
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The result for a random deteetion veetor ean be obtained by replacing (2iV — M) with 1 in the above 
expression, and the corresponding upper bound becomes 

^ ^ 2M[em' + 20w] (r“+2 - 

- (r2-r2)(a + 2) 

which is exactly the same result as the one shown in Lemma [H except for the extra term M which was 
introduced by upper bounding the outage probability in (|9^ . Hence, the proof is completed. 

References 

[1] Y. Saito, Y. Kishiyama, A. Benjebbour, T. Nakamura, A. Li, and K. Higuchi, “Non-orthogonal multiple access (NOMA) for cellular 
future radio access,” in Proc. IEEE Vehicular Technology Conference, Dresden, Germany, Jun. 2013. 

[2] Y. Saito, A. Benjebbour, Y. Kishiyama, and T. Nakamura, “System level performance evaluation of downlink non-orthogonal multiple 
access (NOMA),” in Proc. IEEE Annual Symposium on Personal, huloor and Mobile Radio Communications (PIMRC), London, UK, 
Sept. 2013. 

[3] Z. Ding, Z. Yang, R Fan, and H. V. Poor, “On the performance of non-orthogonal multiple access in 5G systems with randomly 
deployed users,” IEEE Signal Process. Letters, vol. 21, no. 12, pp. 1501-1505, Dec 2014. 

[4] 3rd Generation Partnership Project(3GPP), “Study on downlink multiuser superposition transmission for LTE,” Mar. 2015. 

[5] “5G radio access: requirements, concepts and technologies,” NTT DOCOMO, Inc., Tokyo, Japan, 5G Whitepaper, Jul. 2014. 

[6] “Proposed solutions for new radio access,” Mobile and wireless communications Enablers for the Twenty-twenty Information Society 
(METIS), Deliverable D.2.4, Feb. 2015. 

[7] S. Timotheou and 1. Krikidis, “Fairness for non-orthogonal multiple access in 5G systems,” IEEE Signal Process. Letters, vol. 22, 
no. 10, pp. I647-I65I, Oct. 2015. 

[8] V. Kalokidou, O. Johnson, and R. Piechocki, “A hybrid TIM-NOMA scheme for the SISO broadcast channel,” in Proc. IEEE 
International Conference on Commun. (ICC), Jun. 2015. 

[9] M. Al-Imari, P. Xiao, M. A. Imran, and R. Tafazolli, “Uplink non-orthogonal multiple access for 5G wireless networks,” in Proc. Ilth 
International Symposium on Wireless Communications Systems (ISWCS), Barcelona, Spain, Aug 2014, pp. 781-785. 

[10] H. Marshoud, V. M. Kapinas, G. K. Karagiannidis, and S. Muhaidat, “Non-orthogonal multiple access for visible light communications,” 
arXiv preprint. Available on-line at arxiv.org/abs/1504.00934. 

[11] J. Choi, “Minimum power multicast beamforming with superposition coding for multiresolution broadcast and application to NOMA 
systems,” IEEE Trans. Commun., vol. 63, no. 3, pp. 791-800, Mar. 2015. 

[12] Q. Sun, S. Han, C.-L. I, and Z. Pan, “On the ergodic capacity of MIMO NOMA systems,” IEEE Wireless Commun. Letters, (to appear 
in 2015). 

[13] Z. Ding, F. Adachi, and H. V. Poor, “The application of MIMO to non-orthogonal multiple access,” IEEE Trans. Wireless Commun., 
(to appear in 2015) Available on-line at arXiv: 1503.05367. 

[14] K. Higuchi and Y. Kishiyama, “Non-orthogonal access with random beamforming and intra-beam SIC for cellular MIMO downlink,” 
in Proc. IEEE Vehicular Technology Conference, Las Vegas, NV, US, Sept. 2013, pp. 1-5. 

[15] N. Lee, J.-B. Lim, and J. Chun, “Degrees of the freedom of the MIMO Y channel: Signal space alignment for network coding,” IEEE 
Trans Inform. Theory, vol. 56, pp. 3332 - 3343, Jul. 2010. 

[16] Z. Ding and H. Poor, “A general framework of precoding design for multiple two-way relaying communications,” IEEE Trans. Signal 
Process., vol. 61, no. 6, pp. 1531-1535, Mar. 2013. 

[17] M. Peng, C. Wang, J. Li, H. Xiang, and V. Lau, “Recent advances in underlay heterogeneous networks: Interference control, resource 
allocation, and self-organization,” IEEE Commun. Surveys Tuts., vol. 17, no. 2, pp. 700-729, Second quarter 2015. 

[18] “C-RAN: The road towards green RAN,” China Mobile Res. Inst., Beijing, China, Oct. 2011, White Paper, ver. 2.5. 

[19] Z. Ding, P. Fan, and H. V. Poor, “Impact of user pairing on 5G non-orthogonal multiple access,” IEEE Trans. Vehicular Technology, 
(submitted) Available on-line at arXiv: 1412.2799. 

[20] M. Haenggi, Stochastic Geometry for Wireless Networks. Cambridge University Press, Cambridge, UK, 2012. 

[21] J. Venkataraman, M. Haenggi, and O. Collins, “Shot noise models for outage and throughput analyses in wireless ad hoc networks,” 
in Proc. IEEE Military Communications Conference, Washington, DC, USA, Oct. 2006. 

[22] N. Lee and J.-B. Lim, “A novel signaling for communication on MIMO Y channel: Signal space alignment for network coding,” in 
Proc. IEEE International Symposium on Inform. Theory (ISIT-09), Jul. 2009. 

[23] Z. Ding, T. Wang, M. Peng, and W. Wang, “On the design of network coding for multiple two-way relaying channels,” IEEE Trans. 
Wireless Commun., vol. 10, no. 6, pp. 1820-1832, Jun. 2011. 

[24] M. Rupp, C. Mecklenbrauker, and G. Gritsch, “High diversity with simple space time block codes and linear receivers,” Proc. IEEE 
Global Commun. Conf, vol. 2, pp. 302-306, San Francisco, USA, Dec. 2003. 

[25] C.-H. Liu and J. Andrews, “Multicast outage probability and transmission capacity of multihop wireless networks,” IEEE Trans. Inform. 
Theory, vol. 57, no. 7, pp. 4344^358, Jul. 2011. 

[26] 1. Krikidis, “Simultaneous information and energy transfer in large-scale networks with/without relaying,” IEEE Trans. Commun., 
vol. 62, no. 3, pp. 900-912, Mar. 2014. 

[27] 1. S. Gradshteyn and 1. M. Ryzhik, Table of Integrals, Series and Products, 6th ed. New York: Academic Press, 2000. 




